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Abstract 

The gravitino is a promising super symmetric dark matter candidate which does 
not require exact i?-parity conservation. In fact, even with some small i?-parity 
breaking, gravitinos are sufficiently long-lived to constitute the dark matter of the 
Universe, while yielding a cosmological scenario consistent with primordial nucle- 
osynthesis and the high reheating temperature required for thermal leptogenesis. 
In this paper, we compute the neutrino flux from direct gravitino decay and gauge 
boson fragmentation in a simple scenario with bilinear i?-parity breaking. Our 
choice of parameters is motivated by a proposed interpretation of anomalies in the 
extragalactic gamma-ray spectrum and the positron fraction in terms of gravitino 
dark matter decay. We find that the generated neutrino flux is compatible with 
present measurements. We also discuss the possibility of detecting these neutri- 
nos in present and future experiments and conclude that it is a challenging task. 
However, if detected, this distinctive signal might bring significant support to the 
scenario of gravitinos as decaying dark matter. 



1 Introduction 



The question of the nature of dark matter is still one of the unsolved mysteries in mod- 
ern cosmology. Many particle candidates have been put forward, but until now the only 
dark matter evidence we have is based on the gravitational interaction. In the context 
of supersymmetry with conserved i?-parity, one naturally encounters one of the most 
favoured solutions to the dark matter problem. There, the lightest supersymmetric par- 
ticle (LSP) is stable and can be a successful dark matter candidate if it is neutral and 
weakly interacting like the neutralino. The neutralino is the most thoroughly studied dark 
matter candidate and will be tested in the near future in accelerator, direct detection 
and indirect detection experiments [T]. 

On the other hand, it is also possible that the dark matter interacts only gravita- 
tionally, and supersymmetry also offers candidates of this type. A prominent example is 
the gravitino, the superpartner of the graviton, which was the first supersymmetric dark 
matter candidate proposed [2]. It is one of the most elusive dark matter candidates due 
to its extremely weak interactions. In fact, as part of the gravity multiplet, all gravitino 
interactions are suppressed either by the Planck scale (for the spin-3/2 component) or 
by the supersymmetry-breaking scale (for the Goldstino component). 

Usually, having an LSP with such extremely weak interactions poses a severe problem 
to Big Bang nucleosynthesis, as it makes the next-to-lightest supersymmetric particle 
(NLSP) so long-lived that it decays during or after the formation of the primordial nuclei, 
typically spoiling the successful predictions of the standard scenario [3] . Moreover, if the 
NLSP is charged, the formation of a bound state with ^He catalyzes the production of 
^Li [1] leading to an overproduction of ^Li by a factor 300-600 [5j. One way to avoid 
these constraints is to lower the scale of supersymmetry breaking, thus enhancing the 
Goldstino interactions. However, the reheating temperature of the Universe has to be 
lowered accordingly in order to avoid overclosure ^ and is then in conflict with the 
minimal value required by thermal leptogenesis in order to explain the baryon asymmetry 
of the Universe [7j. 

It was recently proposed in [8] that these problems are automatically solved if a 
small breaking of i?-parity is introduced in the model. Even though in the presence of 
i?-parity violation the neutralino LSP is too short-lived to play the role of dark matter, 
the gravitino LSP can still have a sufficiently long lifetime, which is typically many 
orders of magnitude greater than the age of the Universe due to the suppression of the 
decay rate by the Planck scale and the small i?-parity violating couplings [9]. In this 
scenario, the NLSP population in the early Universe quickly decays into Standard Model 
particles via i?-parity violating interactions. Apart from the presence of a rather inert 
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population of gravitinos, produced thermally from superpartner scatterings at reheating, 
the cosmology then reduces to the non-supersymmetric case well before the synthesis of 
primordial nuclei. 

An intriguing feature of the scenario with i?-parity breaking is that gravitino dark 
matter is not necessarily invisible anymore, since it will decay into Standard Model 
particles at a very slow rate. Since the huge number of gravitinos in our own Galaxy, 
as well as in nearby galaxies and clusters, may compensate for the highly suppressed 
decay rate, this opens up the possibility of observing the dark matter decay products 
as an anomalous contribution to the diffuse gamma-ray flux [101 [H] or the cosmic-ray 
antimatter fluxes [I2l [13]. Indeed, anomalous excesses have been observed both in the 
diffuse extragalactic gamma-ray spectrum and in the positron fraction in similar energy 
ranges. It has been pointed out that the decay of gravitino dark matter with a lifetime 
of ~ 10^^ s and a mass of ~ 150 GeV can account for both of these excesses at the same 
time fl2[ [13] . This motivates our study of the corresponding neutrino spectrum for the 
same choice of parameters, both as a consistency check and to find out whether in this 
scenario an anomalous contribution to the neutrino flux may be expected in present and 
future neutrino experiments. 

This paper is organised as follows: in the next section we will briefly review bilinear 
i?-parity violating models and discuss the resulting gravitino decay modes. In Section [3] 
we will present the neutrino spectrum from dark matter decay and describe its main 
features. In Section [H we will then give the neutrino flux as a function of the gravitino 
lifetime and mass both for neutrinos from our own Galaxy and from diffuse extragalactic 
sources and consider the effect of neutrino oscillations on the signal expected at the 
Earth. In Section [5] we will discuss the different neutrino backgrounds in the energy 
range we are interested in and compare them to our signal. In Section [6| we will propose 
strategies to disentangle the signal from the background, compare the result to present 
neutrino data from Super-Kamiokande. We will then discuss the feasibility of detection 
in future detectors in Section [7| and conclude in Section [HI 



2 Bilinear i?-parity breaking 

We consider here a model of the type described in detail in [8], where the breaking of 
i?-parity is related to the breaking of B — L. In this class of models, the lepton number 
violation in the superpotential is encoded in the bilinear term |14] 

Wf^ = fXiLiHu . (2.1) 
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It can be shown that the _R-parity violating couphng is suppressed compared to the 
coefficient fi of the i?-parity conserving bilinear term H^Hd by /ij ~ ^v'^^^/Mp, with 
vb-l and Mp being the scale of B — L breaking and the Planck mass, respectively. On 
the other hand, the induced baryon number violation is negligible, the corresponding 
Yukawa coupling being further suppressed by the ratio of the gravitino mass to the 
Planck mass and a higher power of the scale oi B — L breaking over the Planck scale, 
A" ~ ms/2V%_i/Mp. Then, if the scale oi B — L breaking is low enough, the constraints 
on proton stability are satisfied and the lepton number violating interactions are kept 
out of equilibrium before the electroweak phase transition, thus preventing the erasure 
of any previously generated baryon asymmetry. This condition requires approximately 
< 10"^"'' [I5], which translates into vb~l ^ 10^'* GeV, although this could be 
circumvented for some specific fiavour structures in the Yukawa couplings. Since the 
present bounds on the gravitino lifetime are already constraining the parameters at this 
order from observations of both gamma rays [TOl [11] and antiprotons ^121 |T3] , this last 
option is probably not viable. 

Apart from the supersymmetric term above, the corresponding soft bilinear super- 
symmetry-breaking term BifiiLiHu arises in the Lagrangian. Since the Bi and /i, terms 
are not usually aligned at the weak scale, a non-vanishing vacuum expectation value 
(v.e.v.) is generated along the sneutrino field direction explicitly breaking lepton number 
and generating not only one neutrino mass, but also non-vanishing mixings between 
neutralinos and neutrinos, as well as between charginos and charged leptons [16]. Such 
mixings are responsible for the two-body decays of the gravitino into gauge boson and 
neutrino, which are the main source of neutrino fiux in our scenario. These decays are 
also possible at the one-loop level if only trilinear i?-parity breaking terms are considered, 
and even in this case they can dominate in part of the parameter space [17]. Since the 
neutralino-neutrino mixing takes place along the Zino component, the branching ratios 
into the different gauge boson channels are fixed by the neutralino mixing matrix once 
the gravitino mass is specified. 

If the gravitino is lighter than the massive gauge bosons, the dominant decay channel 
is the two-body decay into monoenergetic photon and neutrino. Since all the observed 
Yukawa couplings are largest for the third generation, it is reasonable to assume that the 
i?-parity breaking couplings are also largest for the third generation. Thus, the sneutrino 
acquires a v.e.v. only along the Ur direction, and the gravitino will dominantly decay into 
neutrinos with r fiavour. We will see later that this is not a crucial assumption, since 
neutrino oscillations change any pure neutrino fiavour into a mixed state. In particular, 
due to maximal atmospheric mixing, the fiux of tau and muon neutrinos turns out to be 
identical. 
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If instead the gravitino is heavier than the electroweak gauge bosons, the decay into 
these particles is favoured, since the sneutrino has electroweak charge. In this case, neu- 
trinos are produced not only in the direct decay into Z^v^^ but also in the fragmentation 
of the and bosons and the decay of the leptons, thus adding a continuous 
component to the spectrum. In this paper, we will mostly consider the second case of 
a heavier gravitino, motivated by the interpretation of anomalies in other channels as 
explained above. Furthermore, at lower gravitino masses, the detection of any signal is 
much more difficult due to the lower neutrino yield and the smaller signal-to-background 
ratio. 



2.1 Gravitino Decay 

In the models with bilinear i?-parity breaking and a non-zero sneutrino v.e.v. along the 
Vr direction, the main decay channels for the gravitino are: 



^3/2 
^3/2 
^3/2 
^3/2 



(2.2) 



The first decay is practically always allowed, while the next two are open only for a 
gravitino mass above the threshold for or Z^ production. At even higher masses, 
the decay into Higgs boson and neutrino opens up via the Higgsino component of the 
neutralino and the i?-parity violating Higgs-sneutrino mixing. We will consider here the 
case of a large Higgs mass parameter /i, where the lightest Higgs is Standard Model-like 
and the other Higgses decouple. 

The decay widths for these processes can be computed from the interaction La- 
grangian of a gravitino with a gauge boson and a gaugino or with the two chiral Higgs 
multiplets with the insertion of a sneutrino v.e.v. [18j. The decays into electroweak gauge 
bosons arise both from 3-vertices and from the non-abelian 4-vertex, so their structure is 
more complicated than for the abelian sector. The results depend on the gaugino mixing 
matrices and are given by the following: 
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(2.3) 
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with = {i^t)/v, where v = 174 GeV is the Higgs v.e.v.. Assuming ^ 1, the mixing 
matrix elements are given by 



a=l 
^ C9„ C% 
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U-z-z-^MzJ:^^, (2.4) 
L Cw+,.+ c~_. _ + h.c. 



J J Mm/ ^w+xf^w 



'Xa 



2 ^ M± 

Here, Cij are the elements of the unitary matrices that diagonalise the neutralino/chargino 
mass matrices. The kinematical factors are given by 
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The non-abelian 4-vertex computed here was previously neglected in [TT]. Note that 
our results for the decay rates do not agree exactly with those given in |T3]: for the 
interference terms proportional to ji we find a negative sign and a larger coefficient. As 
a result, the branching ratios into the massive gauge boson channels are slightly smaller 
than in [13], while the branching ratio into photon-neutrino is larger. 

The decay into the lightest Higgs, on the other hand, is given by 



^>3/2 .... ..... . . 



3847rMp 



i2 



2 2 

m~ — mf^ 



2 



(2.6) 



in the limit where the hghtest Higgs is given by h = (Re[Hy] sin (3 + Re[Hd] cos/3), 
with 

U^rz^Mzi:"-^^^, (2.7) 

and expressing the soft supersymmetry-breaking Higgs-sneutrino mixed mass term 
through the sneutrino v.e.v. as m|^^T-/v^- 

Note that since the decay into the Higgs is strongly suppressed by the phase space 
factor and since we will consider in the following mostly a gravitino mass of 150 GeV, 
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while m/j > 114 GeV for a Standard Model-like Higgs, this channel is negligible in our 
case, as can be seen from Fig. [H We will therefore ignore it in the following discussion, 
although it will be included in the numerical results and in the figures. 



3 Neutrino Energy Spectrum 



The injection spectrum of neutrinos from gravitino decay is composed of a series of 
contributions. Firstly, the two-body gravitino decay into a photon and a tau neutrino 
produces a monoenergetic line at half the gravitino mass: 



dN,. 



- (V^3/2 ^ l^r) ^5(E 



"^3/2 j 



dE ^'-"^ - '/ V 2 
Additionally, the decay into Z'^Ur produces a second line at an energy 
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which is not completely monoenergetic due to the width of the Z° boson. Instead, it is 
described by a normalised Breit-Wigner profile: 



-1 
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where 



VtZ 



DM. 



m3/2 



(3.3) 



(3.4) 



Furthermore, the decay into hur produces a similar line, with a differential energy 
spectrum denoted by dNu^h/dE. This line is less prominent due to the suppressed branch- 
ing ratio into the Higgs decay channel. Lastly, the fragmentation of the massive gauge 
bosons produces a continuous spectrum of neutrinos in all flavours. We have simulated 
the fragmentation of the gauge bosons with the event generator PYTHIA 6.4 [19j and 
extracted the spectra in the different neutrino flavours for the W^, Z° and h channels, 
which we denote by dN^ /dE, dN^^/dE and dN^^^^/dE, respectively. The leptonic decays 
W —>■ lu, —>■ vv and h ^ Iv also produce, in the rest frame of the decaying particle, 
monoenergetic neutrinos in all flavours. However, due to the boost of the gauge bosons 
in different directions, the lines smear out almost completely in the Earth's rest frame, 
giving just an additional contribution to the continuous part of the spectrum. Taking the 
various decay channels into account, the total spectra for the different neutrino flavours 
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Figure 1: Branching ratios of the different gravitino decay channels as a function of 
the gravitino mass for the case of large n parameter (here taken as /i = 10 TeV) and 
decoupling of the heavy Higgses. The light Higgs mass is taken to be 115 GeV. The 
gaugino mass parameters are taken so as to satisfy the unification relation and to give 
Ml = 1.5 7713/2 at the electroweak scale. The sneutrino mass parameter is fixed at 2m3/2 
and tan f3 = 10. We see that above the mass threshold, the dominant line from Z^Vt 
always has a branching ratio larger than 25 %. 



are given by: 



dE ' ' dE ' ' dE ' ' dE ' 

-dE = ^^^^'^^^ - + ^^(^ -dE- + ^^(^'^^^ -dE- 

, + dN^ , „ , dN^ , ^ dN^. 



(3.5) 



' ^ ciE ^ ' dE 

The branching ratios for the different decay channels can be straightforwardly com- 
puted from the decay widths, Eq. fl2.3p and (12.61) . They turn out to depend mainly on the 
gravitino mass, with a milder dependence on the ratio between the gaugino masses M\ 
and M2 at the electroweak scale. For illustration, we show in Fig. [T]the branching ratios 
as a function of the gravitino mass for the case of a large [i parameter, unified gaugino 
masses with Mi = 1.5m3/2 and = 2 7723/2 at the electroweak scale, and tan/5 = 10. 
In addition, we list in Tab. [H the branching ratios into the different decay channels for a 
number of specific gravitino masses. Note also that in all of the parameter space above 
the Z^ threshold, at least two neutrino lines are present with more than 1 % branching 
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BR(7z/.) 




BR(ZV^) 




lOGeV 


1 








85GeV 


0.53 


0.47 






100 GeV 


0.08 


0.83 


0.09 




150 GeV 


0.01 


0.70 


0.28 


0.01 


250 GeV 


0.003 


0.60 


0.29 


0.11 


1000 GeV 


0.0002 


0.51 


0.25 


0.24 



Table 1: Branching ratios into the different gravitino decay channels for a number of 
specific gravitino masses. The various parameters are chosen as in Fig. [H 
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Figure 2: Neutrino injection energy spectra from gravitino decay for the different flavours 
in the case of ^3/2 = 150 GeV. The line is clearly visible at an energy of 47 GeV. We 
consider here a Higgs mass of 115 GeV, giving a peak at ~ 30 GeV. 



ratio. Actually, sufficiently above the Higgs threshold, the neutrino lines from the Higgs 
and decay channels have comparable strength. 

The energy spectra for the different neutrino flavours, which were obtained as de- 
scribed above, are shown in Fig. [2l In this case, a gravitino mass of ^3/2 = 150 GeV was 
used, and for that value, all three lines mentioned above are visible in the spectrum. We 
see that the tau neutrino spectrum shows a very characteristic signature of two or three 
distinctive peaks (three in the region above the lightest Higgs threshold, as in this case) 
in addition to a continuum which is suppressed at low energies. The spectra of the other 
two flavours are very similar to each other, consisting only of a continuum contribution 
practically following a power law behaviour oc E''^ below the sharp threshold at half the 
gravitino mass. 
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4 Neutrino Fluxes 



In this section we will consider the diffuse neutrino baclcground at the Earth arising from 
the neutrino spectrum from gravitino dark matter decay as discussed above. The diffuse 
neutrino flux has two sources: the decay of gravitinos at cosmological distances and the 
decay of gravitinos in the Milky Way halo. The former contribution is perfectly isotropic, 
while the latter has a mild dependence on the Galactic coordinates [lOl [TT] . 

The decay of gravitinos at a comoving distance xi^)} where z denotes redshift, pro- 
duces a neutrino flux with a redshifted energy spectrum dN^/d{yE), with y = 1 + z. 
Making use of the variation of comoving distance with respect to redshift in a matter- 
and dark energy-dominated Universe, dx/dz = (l + z)~^/^/(ao-ffo\/^A/(l + «:(! + z)~^)), 
with Oo and Hq being the present cosmic scale factor and the Hubble parameter, respec- 
tively, and K = Q/i^/Qm ~ 3 being the ratio between the vacuum and matter density 
parameters, it is straightforward to show that the neutrino flux received at the Earth 
reads: 

oo 

dJ f dN 

Aejdy4^^^^=, (4.1) 



dE "d{yE)./TT^ 
where most neutrinos come from very low redshifts. In this equation, 

(4,2) 

where we have taken the gravitino density to be equal to the cold dark matter den- 
sity, VL^i2h? = 0.1, and the other constants are the critical density pc = 1.05 h"^ x 
10~^ GeVcm^^, the matter density parameter Qm = 0.25 and Ho = 100 hkms~^ Mpc^^ 
with h = 0.73. 

In addition to the extragalactic signal there exists a slightly anisotropic neutrino flux 
stemming from the decay of gravitinos in the Milky Way halo. The energy spectrum is 
given by 

C^-'halo - dN^ 

(4.3) 

where the intensity of the flux, Ahaio, depends on the direction of observation. It is 
proportional to the line-of-sight integration over the halo density profile phaio, being 
defined as 

Aalo = 1 / PhM^)df. (4.4) 



47rr3/2m3/2 j 

l.o.s. 

For our numerical analysis, we will adopt the spherically symmetric Navarro, Frenk and 
White (NEW) profile [20]: 

Phalo(r) = , I I ^12 ' (4-5) 

(r/re)[l + (r/rJJ^ 
10 



with po = 0.26 GeV/cm^ and Vc = 20kpc. The normahsation is chosen such that pir^) = 
0.3GeV/cm^, where Tq = 8.5 kpc is the distance of the Sun from the Galactic centre. 
Our conclusions will turn out to be rather insensitive to the particular choice of the halo 
profile due to the linear dependence of the neutrino fluxes on the dark matter density 
along the line-of-sight and the fact that we will integrate the signal over the whole sky 
excluding the Galactic disk. 

After being produced in gravitino decays, neutrinos propagate while undergoing 
flavour oscillations. Since neutrinos typically travel very long distances before reaching 
us, the conversion probabilities are [21j : 

P(z/e z/^) = ^ (s23 sin^ 26*13 + C23 sin^ 2612) , 

Pii^e ^^r) = \ (4 sin' 2^13 + sl^ siu' 2^12) , (4.6) 
P(z/^ ^ z^r) = ^ (ci3 sin' 26*23 - S23C23 sin' 26*12) , 
while the survival probabilities are 

P(z/e ^ Z/e) = 1 - ^ (sin' 2^13 + 4 siu' 2^12) , 

P{^, ^i^,) = l-\ (44^^3(1 - cl,sl^) + 4 Sin' 2^12) , (4.7) 

P{u^ ^ ly^) = 1 - i (4 44(1 - 44) + ^3 sin^ 2^12) , 

with Sij = sin 6ij and Cij = cos 6ij. Inserting into these equations the experimental best 
fit values for the neutrino mixing angles sin 012 = 0.304, sin' ^23 = 0.5 and sin'^is = 
0.01 [22J, we finally obtain 

P(z/e^z/e) = 0.56, 

P{Ue^Ur) =0.22, (4.8) 

P{Ur ^ Ur) = 0.39 . 

Thus, even when the primary neutrino flux is originally mainly composed of tau neu- 
trinos, the flavour oscillations during the propagation will produce comparable fluxes of 
electron, muon and tau neutrinos due to the large neutrino mixing angles. In particular, 
due to the maximal atmospheric mixing angle, the fluxes of muon and tau neutrinos are 
expected to be essentially identical. 

The fluxes for the different neutrino flavours and their extragalactic and halo con- 
tributions are shown in Fig. [31 In this plot, a band of ±10° around the Galactic disk 
has been removed, and the spectrum is shown with a 10 % energy resolution in order to 
take the finite energy resolution of the detector into account. Note that even with this 
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E (GeV) 

Figure 3: Total neutrino fluxes, as well as the extragalactic and halo contributions for 
the different neutrino flavours after propagation using an energy resolution of 10 %. The 
gravitino mass and lifetime are chosen to be m^i2 = 150 GeV and — 10^^ s. 

optimistic assumption for the energy resolution, the lines from the decay into 7//^ and 
huj. become practically indistinguishable from the continuum, whereas the line from the 
decay into Z^Uj. can be resolved. Its position could allow a determination of the gravitino 
mass, even without determining the endpoint of the spectrum. 

5 Neutrino Backgrounds 

The detection of a possible neutrino signal from gravitino decay is hindered by consid- 
erable neutrino backgrounds. Namely, in the energy range of interest, there exist large 
background neutrino fluxes produced by interactions of cosmic rays with the Earth's 
atmosphere or with the solar corona, as well as neutrino fluxes from distant Galactic 
sources. Let us briefly discuss these backgrounds separately. 

The collision of energetic cosmic rays with nuclei in the upper atmosphere produces 
showers of hadrons, mostly pions, that in turn produce in each decay two muon neutrinos 
and one electron neutrino. Neutrinos arrive from all directions at the detector site after 
propagating a distance ranging between ~ 10 and 12800 km while undergoing flavour 
oscillations. The electron and muon neutrino fluxes have been carefully computed, as- 
suming massless neutrinos, by Battistoni et al. with the Monte Carlo simulation package 
FLUKA [23]. The theoretical uncertainty of the neutrino flux is estimated to be better 
than a 20 % over the energy range of interest for this paper [2l] . 
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The effects of neutrino oscillations on the flavour composition of these fluxes can 
easily be included using the following expression for the conversion probability of muon 
neutrinos into tau neutrinos: 

Piiy>. ^r) = sin2 2^atmsin2 M.27 gj^gv] ) ■ ^^-^^ 

In this expression, E is the neutrino energy and L is their propagation length after being 
produced in the atmosphere, which is given by 

L = v/(i?e cos ey + 2R^h + - cos 6 , (5.2) 

with i?0 = 6371 km being the mean Earth radius and h = lb km the mean altitude 
at which atmospheric muon neutrinos are produced. Moreover, the neutrino parameters 
relevant for the atmospheric oscillations are sin^ 26'atm = 1, I ^"^atml — 2-4 x 10~^ eV^ [22j. 

In addition to the flux of tau neutrinos originating from the conversion of muon 
neutrinos, there exists an intrinsic contribution from the decay of charmed particles 
produced in the atmosphere, coming from all directions, which has a size about 10^ 
times smaller than the flux of electron and muon neutrinos from pion decay. This intrinsic 
contribution has been computed by Pasquali and Reno p5] and can be parametrised as 



log 



10 



^3 dJ,^ I ( GeV^ 



dE / V cm^ s sr 



-A^Bx- Cx^ - L>x^ (5.3) 



where x = log^o (^[GeV]), A = 6.69, B = 1.05, C = 0.150 and D = -0.00820. The 
next-to-leading order QCD calculation also shown in their paper gives lower fluxes for 
energies below several TeV and is therefore less conservative. 

Analogous to the production of neutrinos in the Earth's atmosphere, neutrinos are 
produced in the solar corona by cosmic-ray collisions. This neutrino flux has been studied 
by Ingelman and Thunman in [26j, who found that the flux of electron and muon neu- 
trinos intergrated over the solar disk can be described by the following parametrisation: 

■'^ A^o '^IQ^STl (GeVcn.-r-. (5.4) 



dE 1 + A{E[GeV]) 

which is valid for 10^ GeV < E < 10^ GeV. The numerical values of the coefficients A'^o, 
A and 7 can be found in Tab. [2]for x = 1/^ + i^e, i^ti + i^ti- The electron and muon neutrinos 
and antineutrinos produced in the solar corona oscillate during their propagation to the 
Earth. In view of the long distance travelled, the conversion and survival probabilities can 
be averaged, and the fluxes at the Earth in the different flavours can be straightforwardly 
calculated from Eq. (14.81) . 
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Flavour 


No 


7 


A 


Z/e + Pe 


7.4 X 10"^ 


2.03 


8.5 X 10^^ 




1.3 X 10^5 


1.98 


8.5 X 10-6 



Table 2: Values for the parametrisation of the corona electron and muon neutrino flux. 

Lastly, the fluxes of tau neutrinos that originate from Galactic sources are discussed 
by Athar, Lee and Lin in [27J. For the tau neutrino flux from the Galactic plane in the 
presence of neutrino oscillations they find the parametrisation 



9 X 10-''(GeVcm2ssr)-^ {E[GeV]y^-^\ (5.5) 



dE 

which is valid in the energy range 1 GeV < E < 10^ GeV. 



6 Detection Prospects 

6.1 Comparison with the Signal 

The full-sky signal for the neutrinos from gravitino decay is shown in Fig. H] together with 
the results for the atmospheric background from FLUKA. The signal lies several orders of 
magnitude below the expected atmospheric background for all flavours. Therefore, we find 
that the interpretation of the EGRET and HEAT anomalies in terms of gravitino decay 
is compatible with neutrino flux measurements, as it does not lead to an overproduction 
of neutrinos. 

Going beyond this consistency check, we examine in the following the possibility 
of detecting this exotic contribution in neutrino experiments. Due to the low signal-to- 
background ratio, the signal cannot be detected directly. It will therefore be necessary 
to find strategies for effectively reducing the background in order to have any chance of 
detecting the signal. As is apparent from Fig. HJ the tau neutrino channel appears to be 
the most promising of the three flavours, since it has the lowest background. In general, 
the neutrino spectrum from gravitino decay has some very specific features that could 
allow to distinguish it from the featureless backgrounds, but the question is whether 
neutrino detectors will be able to reach sufficient sensitivity to resolve these features. 

6.2 Electron and Muon Neutrinos 

For the electron and muon neutrinos, which are more easily detected in neutrino ob- 
servatories, the signal-to-background ratio is very low (~ 10~'^~^) for lifetimes that are 
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Figure 4: Full-sky neutrino fluxes expected at the Super-K site for the same gravitino 
mass and lifetime as in Fig. [31 The energy resolution is taken to be 10 %, which is 
insufficient to resolve the different peaks in the spectrum. 



not already excluded by gamma rays [TOl [TT] or antimatter detection [121 [T3j. As can be 
seen from Fig. HI even the peak of the spectrum is three orders of magnitude below the 
background. This makes distinguishing an exotic signal from the background extremely 
difficult. 

Unfortunately, we could not find a suitable strategy to sufficiently reduce this back- 
ground, e.g. by exploiting directionality. In general the atmospheric neutrino baseline 
is too short for all muon neutrinos to oscillate into another flavour at energies of order 
50-100 GeV. It therefore seems hopeless to try to detect the signal without having prior 
knowledge of the position of the peak in the gravitino decay spectrum. In case information 
on the line is available, e.g. from the detection of a monochromatic gamma-ray line by 
the Fermi Gamma- Ray Space Telescope (FGST, formerly named GLAST) [28], then one 
could perhaps envisage strategies to disentangle the signal from the background. How- 
ever, that would probably require a much better knowledge of the atmospheric neutrino 
flux at the relevant energies and a better energy resolution than is presently available. 



6.3 Tau Neutrinos 

For the tau neutrinos, the signal-to-background ratio is more promising, since it lies above 
~ 10~^ at the peak energy. Moreover, most of the background of tau neutrinos from 
atmospheric oscillations and also the other two subdominant sources of tau neutrinos 
can be effectively reduced by exploiting directionality. In fact, the solar corona neutrinos 
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Figure 5: Down-going tau neutrino fluxes expected at the Super-K site with (left) and 
without (right) the contribution from the Galactic disk and the solar corona for the same 
gravitino mass and lifetime as in Fig. [31 

and the Galactic neutrinos mostly come either from the direction of the Sun or from 
the Galactic plane, which could be excluded from the search to reduce the backgrounds. 
Furthermore, in our energy range the tau neutrinos arising from oscillations of the original 
muon atmospheric neutrinos are mostly generated for oscillation lengths of the order of 
the Earth diameter. This means that we expect a very low tau neutrino background if we 
only consider the flux arriving at the detector from above the horizon. In this way, the 
background of tau neutrinos can be reduced by several orders of magnitude. In Fig. [5l 
we show the fluxes for down-going tau neutrinos at the Super-Kamiokande site. We see 
that in this case, the signal can exceed the simulated background from FLUKA, even 
without cutting away the Sun or the Galactic plane. 

6.4 Observability in Super-Kamiokande 

In Cherenkov telescopes, tau neutrinos can only be observed via charged current (CC) 
interactions and are very difficult to disentangle from the other flavours, since the 
Cherenkov signal is not sufficiently distinctive to allow identification on an event-by- 
event basis. 

The Super-Kamiokande collaboration has developed a statistical method to discrimi- 
nate tau neutrinos from the background of other flavours [29] . Using two different strate- 
gies, namely a likelihood analysis and a neural network, they flnd an efficiency of 43.1 % 
and 39.0 %, respectively, to identify tau neutrinos correctly. However, they still misiden- 
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Figure 6: Events in 100 years of observation in Super-K due to CC interactions of 
and Vr from the atmosphere (dash-dotted) and from gravitino decay (sohd), integrating 
the signal from the zenith direction to the angle Q. The gravitino mass and lifetime are 
chosen as in Fig. [31 



tify 3.8 % and 3.4 %, respectively, of the electron and muon background neutrinos as tau 
neutrinos. Due to the large number of electron and muon neutrino events, the sample of 
tau neutrinos is dominated by misidentified neutrinos. The true tau neutrino events can 
therefore only be extracted on a statistical basis using Monte Carlo methods. In the end, 
the data is found to be consistent with the atmospheric tau neutrino flux and neutrino 
oscillations: The full-sky atmospheric tau neutrino signal results in fact in an expected 
78 events in the Super-K I period and 43 events in the Super-K II period [2S1 ED] ■ 

However, this analysis does not exploit the information about the spectral shape 
of the signal apart from setting a threshold for r lepton production, i.e. Ey^ > 111^. + 
m^/2m„ ~ 3.5 GeV, so this kind of data analysis could certainly be improved to search 
for a signal with a peak above the continuum, as in our case. 

Despite the experimental difficulties, it is worthwhile to examine the theoretically 
expected signal in the tau channel. Fig. [6] shows the expected number of tau neutrino 
and antineutrino events per century of observation at Super-Kamiokande within a zenith 
angle integrated from cos to 1. If only down-going neutrinos are selected, the signal from 
gravitino decays lies above the atmospheric background for higher energies. However, the 
fluxes are extremely low and result in only a few events per century, making it practically 
impossible to discriminate them from the other flavours using statistical methods. 

One detector specifically optimised for measuring tau neutrinos above 17 GeV event- 
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Figure 7: Region of the gravitino lifetime where the hne signals from the two-body 
decays into ^v^ and Z°z/^ overcome the atmospheric background either for all neutrino 
flavours or for down-going tau neutrinos only. We see that considering the tau flavour and 
exploiting directionality can improve the detectability by three orders of magnitude for 
large gravitino masses. We use here the branching ratios shown in Fig.[T]for the two lines. 
Note that in the intermediate region the continuum spectrum from fragmentation 
may be used to close the gap, but we do not consider this possibility here. 



by-event is OPERA [31], which is already active in Gran Sasso and will measure tau 
neutrino appearance in a muon neutrino beam produced at CERN. Unfortunately, the 
detector's effective mass is more than a factor 10 smaller than that of Super-K and thus, 
even neglecting the issue of directionality, would be able to observe only one event from 
gravitino decays in more than 1000 years. 

We therefore conclude that present detectors are unable to detect the signal, either 
because they do not have a sufficiently good efficiency for identifying tau neutrinos, or 
because they are too small for the low intensity of our signal, or both. 



7 Future Detectors 



In this section, we will briefly discuss the prospect of detection for the future. In order 
to consider the general parameter space for gravitino decay, we show in Fig. [7] the region 
of gravitino lifetime and mass where the neutrino signal in the and gamma peaks is 
equal to the simulated atmospheric background. For the muon/electron neutrino flux, 
this result is similar to that presented in [23 ? except that we are weighting the channels 
with the gravitino branching ratios and that we have only one neutrino produced in the 
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line instead of two. Real limits from neutrino experiments have been considered in the 
past in [33] . 

We are taking here as in [32] a nominal energy resolution of 0.3 in log^o(-^[GeV]) 
around the peak position and we single out the value of the lifetime for which the peak 
is equal to the background. Note that since the signal is proportional to l/r3/2, requiring 
the peak to be larger than the background by a specific factor only rescales the curves 
by the inverse of this factor. 

We clearly see again that the tau neutrino channel in the down-going direction allows 
to constrain the gravitino lifetime a few orders of magnitude better than the muon or 
electron neutrino fluxes. On the other hand, similar plots for the gamma-ray channel [10] 
are even more sensitive and give bounds at the order of 10^'' s for gravitino masses below 
the and Z° thresholds. For masses above 200 GeV, the tau neutrino channel starts 
to compete in sensitivity with the photon channel, if we neglect for the moment the 
difficulties connected with measuring such a low flux and identifying the neutrino flavour. 
At even higher energies the background flux decreases quickly and therefore the signal-to- 
background ratio improves, but the signal rate then also decreases (as 1/7713/2), making 
detection more difficult. 

7.1 Hyper-Kamiokande 

The prospects for Hyper-Kamiokande can be easily obtained by considering that its mass 
is planned to be a factor of 10 (for the 0.5 megaton project) to 20 (for a 1 megaton case) 
larger than Super-Kamiokande. Assuming that the rest of the detector performances are 
unchanged, we expect to find approximately 20-40 events from our signal per century 
from the upper hemisphere]^ 

This number of events might still be too small to allow for statistical analysis. How- 
ever, we expect most of the events to appear within the peak region or near the threshold 
and therefore, an appropriate energy binning, especially optimised after a signal has al- 
ready been detected in gamma rays, could allow to collect a significant number of events 
above the background in a specific energy bin on a shorter timescale. Still, it is clear 
that a sufficiently good energy resolution is a key requirement for singling out the line 
events, and it remains uncertain how and if the tau statistical discrimination analysis 
can be applied to a sample of such few events. 

^ This number could be larger if the fiducial volume of Hyper-Kamiokande is larger than ^1/2 the 
total volume as it is in Super-Kamiokande. 
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7.2 IceCube and km^ Detectors 

Detectors of km^ dimensions have in principle sufficient size to collect enough events 
to detect the signal within a reasonable time span. Even considering that IceCube is 
actually looking downwards and not at the upper hemisphere, from the horizontal direc- 
tion and the proton cross-section we estimate (9(100) events per year for the completed 
experiment. Of course, the effective area depends on the neutrino energy: taking the 
effective area given in [33] for the opposite direction and assuming most of the signal 
is above 100 GeV, we have instead (9(10) events per year. In general, it would be desir- 
able to lower the energy threshold to reach below 100 GeV in order to cover the energy 
range favoured by the EGRET and HEAT anomalies. The combination of IceCube with 
AMANDA already allows to lower the threshold to 30 GeV. Additionally, plans are being 
considered for adding another, denser subdetector at a deeper location to improve the 
sensitivity to dark matter annihilations [3lj. Such a configuration could probably also 
be useful for investigating the present scenario and, more generally, other decaying dark 
matter candidates. 

However, in the case of Cherenkov detectors, the discrimination of tau neutrinos from 
other neutrino flavours is generally difficult, and for IceCube strategies for tau flavour 
identification have been proposed only for neutrinos well above TeV energies [35]. It 
could therefore be more favourable to improve the energy resolution and exploit the 
muon neutrino final state instead. 

8 Conclusions 

We have examined the neutrino spectrum from the decay of unstable gravitino dark 
matter in a scenario with bilinear i?-parity violation. It has been pointed out in the 
recent literature that the decay of gravitino dark matter particles with a lifetime of 
~ 10^^ s and a mass of ~ 150 GeV into massive gauge bosons may account for the 
anomalies observed in the diffuse extragalactic gamma-ray spectrum as measured by 
EGRET as well as in the positron fraction as measured by HEAt|§ Motivated by this 
observation, we have computed the neutrino spectrum for the same choice of parameters 
as a consistency check of this scenario. We find that this spectrum is compatible with 
results from neutrino experiments. 

We have also examined the detectability of this exotic component of the neutrino 
flux to find an independent way to test this scenario. While the signal in the neutrino 
spectrum with two or more distinct peaks, resulting from two-body gravitino decays into 

^The existence of a positron excess seems to be supported by preliminary results from PAMELA [36j. 
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gauge/Higgs boson and neutrino, is very characteristic, it will be challenging to detect 
these features in neutrino experiments. On one side, present neutrino detectors do not 
achieve a sufficiently high energy resolution to resolve the subdominant peaks, and on the 
other side, the event rate is expected to be so small that the background of atmospheric 
neutrinos overwhelms the signal in all flavours. The most promising signal-to-background 
ratio is found in the tau neutrino flavour, especially when analysing only the flux from 
the upper hemisphere since there the atmospheric tau neutrino flux is vastly reduced. 
However, tau neutrinos are difficult to identify in Cherenkov detectors and probably 
only an cvent-by-event identification procedure could allow the signal to be seen with 
such extremely limited statistics. At present, therefore, it is not possible to detect this 
contribution due to technological limitations. 

The ideal detector for testing the present scenario would be one of megaton mass 
with the ability to identify and measure tau neutrinos event by event. Should such a 
detector ever become available, it could be worthwhile to look for this component of the 
neutrino flux by employing strategics for background reduction such as the ones discussed 
here, especially if the anomalous signatures in the positron fraction and the diffuse extra- 
galactic gamma-ray spectrum are con&med by PAMELA and FGST, respectively. The 
detection of a signal in neutrinos compatible with signals in the other indirect detection 
channels would in fact bring significant support to the scenario of decaying dark matter, 
possibly consisting of gravitinos that are unstable due to bilinear i?-parity violation. 
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